A quantum statistical approach to simplified stock 

markets 



F. Bagarello 

Dipartimento di Metodi e Modelli Matematici, Facolta di Ingegneria, 
Universita di Palermo, I - 90128 Palermo, Italy 
E-mail: bagarell@unipa.it 
home page: www.unipa.it\~ bagarell 



Abstract 

We use standard perturbation techniques originally formulated in quantum (statis- 
tical) mechanics in the analysis of a toy model of a stock market which is given in 
terms of bosonic operators. In particular we discuss the probability of transition 
from a given value of the portfolio of a certain trader to a different one. This com- 
putation can also be carried out using some kind of Feynman graphs adapted to the 
present context. 
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I Introduction and motivations 



In some recent papers, [TJ [2j [3], we have discussed why and how a quantum mechanical 
framework, and in particular operator algebras and the number representation, can be 
used in the analysis of some simplified models of stock markets. These models are just 
prototypes of real stock markets because, among the other simplifications, we are not 
considering financial derivatives. For this reason our interest looks different from that 
widely discussed in j3], even if the general settings appear to be very close (and very close 
also to the framework used in [5]). The main reason for using operator algebras in the 
analysis of these simplified closed stock markets comes from the following considerations: 
in the closed market we have in mind the total amount of cash stays constant. Also, the 
total number of shares does not change with time. Moreover, when a trader r interacts 
with a second trader a, they change money and shares in a discrete fashion: for instance, 
r increments his number of shares of 1 unit while his cash decrements of a certain number 
of monetary units (which is the minimum amount of cash existing in the market: 1 cent 
of dollar, for example), which is exactly the price of the share. Of course, for the trader 
a the situation is just reversed. So we have at least two quantities, the cash and the 
number of shares, which change discontinuously as multiples of two fixed quantities. In 
[TJ [21 E] we also have two other quantities defining our simplified market: the price of the 
share (in the cited papers the traders can exchange just a single kind of shares!) and the 
market supply, i.e. the overall tendency of the market to sell a share. It is clear that also 
the price of the share must change discontinuously, and that's why we have assumed that 
also the market supply is labeled by a discrete quantity. 

Operator algebras and quantum statistical mechanics provide a very natural settings 
for discussing such a system. Indeed they produce a natural way for: (a) describing 
quantities which change with discrete steps; (b) obtaining the differential equations for 
the relevant variables of the system under consideration, the so-called observables of the 
system; (c) finding conserved quantities; (d) compute transition probabilities. 

For these reasons we have suggested in [TJ [2j [3] an operator- valued scheme for the de- 
scription of such a simplified market. Let us see why, neglecting here all the many math- 
ematical complications arising mainly from the fact that our operators are unbounded, 
and limiting our introduction to few important facts in quantum mechanics and second 
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quantization which will be used in the next sections. More details can be found, for 
instance, in P, E] and [5], as well as in [TJ [2J [3]. 

Let 7i be an Hilbert space and B(7i) the set of all the bounded operators on 7i. 
Let S be our (closed) physical system and 21 the set of all the operators, which may be 
unbounded, useful for a complete description of S, which includes the observables of S. 
The description of the time evolution of S is driven by a self-adjoint operator H = 
which is called the hamiltonian of S and which in standard quantum mechanics represents 
the energy of S. In the Heisenberg picture the time evolution of an observable X e 21 is 
given by 

X(t) = e iHt Xe- iHt (1.1) 
or, equivalently, by the solution of the differential equation 

^ = ie im [H,X]e- im = i[H,X{t)], (1.2) 

where [A, B] := AB — BA is the commutator between A and B. The time evolution 
defined in this way is usually a one parameter group of automorphisms of 21. The wave 
function \1/ of S is constant in time. 

In the Scrodinger picture the situation is just reversed: an observable X e 21 does not 
evolve in time (but if it has some explicit dependence on t) while the wave function ^ of 
S satisfies the Scrodinger equation i 9 yp = H^/(t), whose formal solution is, if H does 
not depend on time, \&(t) = e~ tHt ^ . 

In our paper a special role is played by the so called canonical commutation relations 
(CCR): we say that a set of operators {ai, aj,l = 1,2,...,L} satisfy the CCR if the 
following hold: 

[a h al\ = 5 in l, [ai, a n ] = [a], a+] = (1.3) 

for all I, n = 1, 2, . . . , L. Here 1 is the identity operator of B{7i). These operators, which 
are widely analyzed in any textbook in quantum mechanics, see [6] for instance, are those 
which are used to describe L different modes of bosons. From these operators we can 
construct hi = ajai and N = Ylt=i which are both self-adjoint. In particular hi is the 
number operator for the 1-th mode, while N is the number operator of S. 

The Hilbert space of our system is constructed as follows: we introduce the vacuum 
of the theory, that is a vector Lp which is annihilated by all the operators af. aifo — 
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for all / = 1, 2, . . . , L. Then we act on (p with the operators a] and their powers: 

<p m , n2 ,..., nL : = 1 , («ir(4r---(4r>o, (1.4) 

ri/ = 0, 1,2, . . . for all /. These vectors form an orthonormal set and are eigenstates of 
both hi and N: hi(p ni>n2i ... >nL = n^ m>ri2) ... in£ and N<p ni ,na,...,n L = N(p ni>n2i ... jriL , where iV = 
Ya=i u i- Moreover using the CCR we deduce that hi {aiip ni , n2 ,...,n L ) = (ni-l)(aWn u n 2 ,...,n L ) 
and hi (al(p nun2v ,. in ^\ = (ni + l)(aj<p nijn2) ___ jnL ), for all /. For these reasons the following 
interpretation is given: if the L different modes of bosons of S are described by the vector 
V 9 ni,n 2 ,...,n i , this implies that n\ bosons are in the first mode, ri2 in the second mode, and 
so on. The operator hi acts on (f nit n 2 ,...,n L an d returns ni, which is exactly the number of 
bosons in the 1-th mode. The operator N counts the total number of bosons. Moreover, 
the operator a>i destroys a boson in the 1-th mode, while a\ creates a boson in the same 
mode. This is why a/ and a\ are usually called the annihilation and the creation operators. 

The Hilbert space 7i is obtained by taking the closure of the linear span of all the 
vectors in ( 11. 4p . 

An operator Z e 21 is a constant of motion if it commutes with H. Indeed in this case 
equation ( 11.21) implies that Z(t) = 0, so that Z(t) = Z for all t. 

The vector y2 ni> n 2 ,...,n L i n d 1-41) defines a vector (or number) state over the algebra 21 

as 

where < , > is the scalar product in Ti. As we have discussed in [U |2], these states may 
be used to project from quantum to classical dynamics and to fix the initial conditions of 
the market. 

The paper is organized as follows: 

In Section II we introduce a new model, slightly different from the one proposed in 
[2] , and we deduce some of its features and the related equations of motion, working in 
the Heisenberg picture. One of the main improvements with respect to [2] is that several 
kind of shares (and not just one!) will be considered here. 

In Section III we adopt a different point of view, using the Scrodinger picture to deduce 
the transition probability from a given initial situation to a final state, corresponding to 
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two different values of the portfolios of the various traders. Since the reader might not be 
familiar with the tools adopted, we will be rather explicit in the derivation of our results. 
Section IV is devoted to the conclusions and to our plans for the future. 

II The model and first considerations 

Let us consider N different traders n, T2, . . ., Tjy, exchanging L different kind of shares a\, 
(T2, . . ., &l- Each trader has a starting amount of cash, which is used during the trading 
procedure: the cash of the trader who sells a share increases while the cash of the trader 
who buys that share consequently decreases. The absolute value of these variations is 
the price of the share at the time in which the transaction takes place. Following our 
previous results we start introducing a set of bosonic operators which are listed, together 
with their economical meaning, in the following table. We adopt here latin indexes to 
label the traders and greek indexes for the shares: j — 1, 2, . . . , N and a — 1, 2, . . . , L. 





the operator and.. 


...its economical meaning 






annihilates a share a a in the portfolio of re- 




a j,a 


creates a share a a in the portfolio of re- 






counts the number of share o~ a in the portfolio of r, 




C i 


annihilates a monetary unit in the portfolio of re- 






creates a monetary unit in the portfolio of re- 




kj = c j c j 


counts the number of monetary units in the portfolio of Tj 




Pa 


lowers the price of the share o~ a of one unit of cash 




pi 


increases the price of the share a a of one unit of cash 




Pa PaPoi 


gives the value of the share a a 



Table 1.- List of operators and of their economical meaning. 



These operators are bosonic in the sense that they satisfy the following commutation 
rules 

[Cj,c[] = 16 jtk , [p a ,p\] = l^a,/? K>, <4 )/9 ] = ISjfiS^p, (2.1) 
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while all the other commutators are zero. 

As discussed in the Introduction, in [21 [3] we have also introduced another set of 
operators related to the market supply which was used to deduce the dynamics of the 
price of the single kind of share considered there. However, the mechanism proposed in 
those papers, thought being reasonable, is too naive and gives no insight on the nature 
of the market itself. For this reason in [51 [H] we have also considered a different point 
of view, leaving open the problem of finding the dynamics of the price and focusing the 
attention on the time evolution of the portfolio of a fixed trader. This is the same point of 
view which we briefly consider in this section, while we will comment on other possibilities 
in the rest of the paper. More precisely, as we have discussed in the Introduction, the 
dynamical behavior of our market is driven by a certain hamiltonian H. We assume that 
H can be written as H — H + H pr i ces , where 

H — H + X Hi, with 

H = J2j,a ™j,a + J2j k j (2.2) 

Hi = Eyo (^iHachf + h.c)j . 

Here h.c. stands for hermitian conjugate, cf a and cj " are defined as in [2], and io^ a , tvj 
and p\ °j are positive real numbers. In particular these last coefficients assume different 
values depending on the possibility of T{ to interact with Tj and exchanging a share a a : 
for instance p^\ = if there is no way for r 2 and r 5 to exchange a share <j\. Notice that 
this does not exclude that, for instance, they could exchange a share 02, so that ^5 7^ 0. 
With this in mind it is natural to put p^ = and pf*- = p^ . 

Going back to (12.21) . we observe that Hq is nothing but the standard free hamiltonian 
which is used for many-body systems like the ones we are considering here (where the 
bodies are nothing but the traders, and the cash). More interesting is the meaning of the 
interaction hamiltonian Hi. To understand Hi we consider its action on a vector like 

t m,i f «at,l ffci f k N |Pi f p l 
°i,i ' ' ' a N,L c i ' ' ' c n Pi ' ' 'Pl /„ x 

f{nj,c,};{kj};{P a } ■= 1 j F=j , , , j <A), (2.3) 

y/nul ■ ■■n N)L \ki\ ■ ■■k L \P 1 \ ■ ■ ■ P L \ 
where <^ is the vacuum of all the annihilation operators involved here, see Section I. 
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Because of the CCR we deduce that the action of a single contribution oiHj, a\ a aj )a cf a c^ *, 
on firij a };{kj};{p a } is proportional to another vector </?{„'. a } ; {^} ; {P^} with just 4 different 
quantum numbers. In particular ny )Q ,, n iiQ , fcj and fc* are replaced respectively by rij >a — 1, 
rij a + 1, kj + P a and fcj — P a (if this is larger or equal than zero, otherwise the vector is 
annihilated). This means that Tj is selling a share o~ a to and earning money from this 
operation. For this reason it is convenient to introduce the following selling and buying 
operators: 

%j,a := a j,aCj j ^j- )Q := a j,a C j " (2-4) 

With these definitions and using the properties of the coefficients pj^ we can rewrite H 



i 



as 



H i = 2 Yl Ptj x U x i,* ^ H = Y1 u i> a + J2 u i" k i + 2X Y1 PiJ x U (2-5) 

i,j,a j,a j i,j,a 

In [5] we have discussed the role of H price which should be used to deduce the time evolution 
of the operators P a , a = 1, 2, . . . , L, and which will not be fixed in these notes. This will 
be justified below, after getting the differential equations of motion for our system. Again 
in [9] we have also shown that H corresponds to a closed market where the money and 
the total number of shares of each type are conserved. Indeed, calling N a := YliLi ^i,a 
and K := Y^iLi^i we have seen that, for all a, [H,N a ] = [H, K] = 0. Hence N a and 
K are integrals of motion, as expected. Of course, something different may happen in 
realistic markets. For instance, a given company could decide to issue more stocks in the 
market or to split existing stocks. The related iVa, say N ao , is no longer a constant of 
motion. Hence H should be modified in order to have [H, N ao ] ^ 0. For that it is enough 
to add in H some source or sink contribution, trying to preserve the self-adjointness of H. 
However, this is not always the more natural choice. For instance, if we want to include 
in the model interactions with the environment (economical, political, social inputs), it 
could be more convenient to use non-hermitean operators, like the generators appearing 
in the analysis of quantum dynamical semigroups used to describe open systems, [TO] . 
This aspect will not be considered here. 
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Let us now define the portfolio operator of the trader 77 as 



ft,(*) = ^2P a (t)hi ta (t) + k l (t). 



(2.6) 



a=l 



This is a natural definition, since it is just the sum of the cash and of the total value of 
the shares that 77 possesses at time t. Once again, we stress that in our simplified model 
there is no room for the financial derivatives. 

Using (II. 2p and the commutation rules assumed so far we derive the following system 
of equations: 



dt 

dkiif) _ 
dt 

dt 



(2.7) 



which, together with their adjoints, produce a closed system of differential equation. 
Notice that these equations imply that ^ =1 P„(t) 4;ni t0 ,(t) + ;§^(A) = ^- This system 
is now replaced by a semi-classical approximation which is obtained replacing the time 
dependent operators P a (t) with L classical fields P a (t) which are deduced by empirical 
data. This is the reason why we were not interested in fixing H pT i ces in H: P a (t) is replaced 
by P a {t) which are no-longer internal degrees of freedom of the model but, rather than 
this, simple external classical fields. Neglecting all the details, which can be found in [9], 
we observe that the first non trivial contribution in A is 



ni >a (t) ■= ^{n jtCl }-,{k j y,{p a }(ni,a( t )) 



= n lja - 8A 2 £ 

ki(t) := aJ {nja y. {kj y. {Pa} (ki(t)) 

" («) 



p$Y M jtl . a X(G ( ?l a (t)) =: n ha + 6n ha (t), 



(2.8) 



= h + 8A 2 e;Li e: = i k • 

where we have defined: 



M^{Qfljfy) =: h + 5h{t), 



M jt i. a = M u . a - Mij. ia , 



Mu, a 



(k J+ p a y. (k l+ p a y. 



- n i)Q ,(l + n h 



(2.9) 



(k t -P a )V 
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and 



Jo ©SU*'K 



Jo e 



Ceo di 

ft 



-ies 

-W,) f*P a (t')dt' 



[2.10) 



10J 



j. a 



The time dependence of the portfolio can now be written as 



n<(*) - ^MkMP^iMt)) = nj(o) + su^t), 



(2.11; 



with 



= ^ n,, a (P a (t) - P a (0)) + p a(t) Sn ha (t) + <JA;,(t). (2.12) 

a=l a=l 

It should be emphasized that, even under the approximations we are considering here, we 
still find J2a=i Fait) ni, a (t) + h(t) = 0. However, we are loosing the conservation laws 
we have discussed before: both TV; and K do not commute with the effective hamiltonian 
which produces the semi-classical version of (12.71) anymore, and therefore they are not 
constant in time. This imposes some strict constraint on the validity of our expansion, 
as we have discussed already in [3] and suggests the different approach to the problem 
which we will discuss in the next section. Again we refer to [9] for more comments on 
these results. 



Ill A time dependent point of view 

In this section we will consider a slightly different point of view. Our hamiltonian has no 
H price contribution at all, since the price operators P a , a — 1, . . . , L, are now replaced from 
the very beginning by external classical fields P a {t), whose time dependence describes, as 
an input of the model, the variation of the prices of the shares. Incidentally, this implies 
that possible fast changes of the prices are automatically included in the model through 
the analytic expressions of the functions P a (t). Hence the interaction hamiltonian Hi in 
(12.51) turns out to be a time dependent operator, Hj(t). More in details, the hamiltonian 



9 



H = H + XHj(t) of the model looks like the one in (12.51) but with the following time- 
dependent selling and buying operators: 

x j>a (t) ■= H«cf ait \ x] >a (t) := at Q c/«W (3.1) 

The L price functions P a {t) will be taken piecewise constant, since the price of a share 
changes discontinuously: it has a certain value before the transaction and (in general) 
a different value after the transaction. This new value does not change until the next 
transaction takes place. More in details, we introduce a time step h which we call the 
time of transaction, and we divide the interval [0,t[ in subintervals of duration h: [0,t[= 
[t ,t 1 [U[t ,t 1 [U[t 1 ,t 2 [- ■ ■ [£m-i,£m[, where t = 0, t x = h, . . ., t M -i = (M - l)h = t - h, 
tM = Mh = t. Hence h = t/M. As for the prices, for a = 1, ... L we put 



Pa,0, t G [to, tl[ 

Pa.li t £ Fl, ^[ 



(3.2) 

, 

Pa,M-l, t G [tM-l,tM[ 



An orthonormal basis in the Hilbert space of the model 7i is now the set of vectors defined 
as 

°1,1 ' ' ' a N,L C l ' ' ' C N /Q q n 

A/Tin ! • ■■n NjL \k 1 \ ■■■k L \ 

where ifo is the vacuum of all the annihilation operators involved here. They differ from 
the ones in (12.31) since the price operators disappear, of course. To simplify the notation 
we introduce a set T = {{nj a }; {%}} so that the vectors of the basis will be simply 
written as tpjr. 

The main problem we want to discuss here is the following: suppose that at t = the 
market is described by a vector ip^ . This means that, since JF = {{n° Q }, {k°}}, at t = 
the trader T\ has n\ x shares of <Ti, n° 2 shares of <T2, . . ., and k° units of cash. Analogously, 
the trader r 2 has n° 2X shares of <7i, n° 22 shares of cr 2 , . . ., and k 2 units of cash. And so 
on. We want to compute the probability that at time t the market has moved to the 
configuration Tf = {{rij a }, {kj}}. This means that, for example, T\ has now n xl shares 
of o"i, n{ 2 shares of 02, . . ., and k{ units of cash. 
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Similar problems are very well known in ordinary quantum mechanics: we need to 
compute a probability transition from the original state </?jr to a final state (fj?., and 
therefore we will use here the standard time-dependent perturbation scheme for which 
we refer to [11]. The main difference with respect to what we have done in the previous 
section is the use of the Schrodinger rather than the Heisenberg picture. Hence the market 
is described by a time-dependent wave function which, for t = 0, reduces to (pjr : 
\l/(0) = tpj? . The transition probability we are looking for is 

P^ f (t):=\<<p^(t)>\ 2 (3.4) 

The computation of Pj? -+j?.(t) is a standard exercise, [TT] . In order to make the paper 
accessible also to those people who are not familiar with quantum mechanics, we give 
here the main steps of its derivation. 

Since the set of the vectors y?jr is an orthonormal basis in Ti, the wave function \l/(t) 
can be written as 

*(t) = X>(t) e-a^V, (3.5) 

T 

where Ejr is the eigenvalue of H defined as 

H Lfjr = EpLfjr, =>• Ejr = ^ L0j, a n j>a + ^ LOjkj. (3.6) 

3,ot j 

This is a consequence of the fact that iff in (13.31) is an eigenstate of H in (12. 2p . Using the 
quantum mechanical terminology, we sometimes call the free energy of (f^. Putting 
(13.51) in (13 .4p . and recalling that < ip^, cpg >= 5^,g, we have 

P^ f (t) := |c^(t)| 2 (3.7) 

The answer to our original question is therefore given if we are able to compute cjrJt) in 
(13.51) . Due to the analytic form of our hamiltonian, this cannot be done exactly. However, 
several possible perturbation schemes exist in the literature. We will adopt here a simple 
perturbation expansion in the interaction parameter A appearing in the hamiltonian ( 12.51 ) . 
In other words, we look for the coefficients in (13. 5p having the form 

Cr (t) = cf (t) + Ac? (t) + \ 2 4 ] (*) + ••• (3.8) 
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Each Cj/(t) satisfies a differential equation which can be deduced as follows: first we recall 
that ty(t) satisfies the Schrodinger equation i 8 ^ — H(i)^(t). Replacing (13.51) in this 



equation and using the orthonormality of the vectors <f F s, we find that 



(3.9) 



Replacing now (13.81) in (13.91) we find the following infinite set of differential equations, 
which we can solve, in principle, up to the desired order in A: 



(3.10) 



The first equation, together with the assumed initial condition, gives cp,{t) = c^(0) = 
$T',Fo- When we replace this solution in the differential equation for c F (t) we get, recalling 
again that ^(0) = <p Fo , 

c «(t) = -i f e^-E^t, < ^ /)ifj ( tl )^ o > dtl (3.H) 
Jo 



Using this in (13. lUp we further get 



cf,{t) 



I 2 e^~ E ^ h^ih) dt^ e^'~ E ^ h^A^) dh, (3.12) 



where we have introduced the shorthand notation 



(3.13) 



III.l First order corrections 

We continue our analysis computing Pj? _>jr.(t) in (13.71) up to the first order corrections 
in A and assuming that J 7 / is different from JF . Hence we have 

2 



(3.14) 
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Using (13.21) and introducing 5E = Ej? f — Ejr Q , after some algebra we get 



5 Eh/2 
6E/2 



M-l 
k=0 



(3.15) 



The computation of the matrix elements hjr f ^ {tk) is easily performed. Indeed, because 
of some standard properties of the bosonic operators, we find that 



t P a ,k to _ r( fc ) 

a i,a a j,a Cj C -T Q ft <y9jr — 1 j • </? (;, 



where 



r<*> .- M + P °*> 1 k ° [ n o (1 + n o) (3 16) 

and J r Qk' a ' > differs from JF only for the following replacements: n° a — > n° a — 1, n° a — > 
n° a +l, k° — > k° + P at k, k° — » k° — P a ^. Notice that in our computations we are implicitly 
assuming that k° > for all i, k and a. This is because otherwise the trader Tj would 
have not enough money to buy a share <r Q ! 
We find that 

hr f ,r (h) = 2][>S ) rg Q < ^,,^,«) > (3.17) 

Of course, due to the orthogonality of the vectors <pp8, the scalar product < cp^ f , <£>_.(»,.*,<*) > 
is different from zero (and equal to one) if and only if n^ a = n° a — 1, n{ a = n° i a + 1, 
H = K ~ Pa,k and kf = k° + P a ,k, and all the other new and old quantum numbers 
coincide. 

For concreteness sake we now consider two simple situations: in the first example 
below we just assume that the prices of the various shares do not change with t. In the 
second example we consider the case in which only few changes occur, and we take M = 3. 

Example 1: constant prices 

Let us assume that, for all k and for all a, P a ^ = P a {tk) = Pa- This means that 
H*/. , Fqu' 01 ^ and the related vectors </?_.(i,j» do not depend on k. Hence f ^ (tk) is also 
independent of k. After few computation we get 



, /'sin(5£t/2)\ 2 ,, , n .\2 



P^ f (t)=y( ^ ' - ) \h^ (0)\ (3.18) 
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to which corresponds a transition probability per unit of time 

Vt^t s = lim-P^^) = 2n\ 2 5(Er f - E Fo ) |/^ (0)| 2 , (3.19) 

t,oo t 

which shows that, in this limit, a transition between two states is possible only if the two 
states have the same free energy. The presence of hj^,jr o (0) in the final result shows, using 
our previous remark, that at the order we are considering here a transition is possible only 
if (yjjr does not differ from cpj? for more than one share in two of the rij^'s and for more 
than P a in two of the k/s. All the other transitions are forbidden. 

Example 2: few changes in the price 

Let us now fix M = 3. Formula (I3.15P can be rewritten as 

JW,M = * I ErS'ffH < n,,^ > + 

+r« Q < tp^ip^) > e ihSE + rg ;Q < ip^tp^ > e 2ihSE )\ 2 (3.20) 

The meaning of this formula is not very different from the one discussed in the previous 
example: if we restrict to this order of approximation, the only possibilities for a transition 
JF — > Tj to occur are those already discussed in Example 1 above. We will see that, in 
order to get something different, we need to go to higher orders in A. In other words, 
even if the prices depend on time, not new relevant features appear in the transition 
probabilities. 

As for the validity of the approximation, let us consider the easiest situation: we have 
constant prices (Example 1) and, moreover, in the summation in (13.171) only one contri- 
bution survives, the one with i , j and «o- Then we have h^ f -f o (to) = 2p-°^r ioJo;Q , < 
Lfjr y_(i j ,ex ) >• Because of (13.181) . and since our approximation becomes meaningless if 

Pjr _>j: (t) exceeds one, it is necessary to have small A, small pf^ and large 8E (if this is 
possible). However, due to the analytic expression for T io ^ . ao , see (13.161) . we must pay at- 
tention to the values of the n^'s and of the k/s, since, if they are large, the approximation 
may likely break down very soon in t. 

Let us now finally see what can be said about the portfolio of the trader n. Since 
we know the initial state of the system, then we know the value of its portfolio at time 
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zero: extending: our orig inal definition in Q233) we have IL,(0) = J2a=i P «(°) H a {0) + k{0) . 
As a matter of fact, the knowledge of <^ implies that we know the time zero value of 
the portfolios of all the traders, clearly! Formula (13.151) gives the transition probability 
from iffr to tpFj. This probability is just a single contribution in the computation of the 
transition probability from a given 11; (0) to a certain Ilj(i), since the same value of the 
portfolio can be recovered at time t for very many different states <p^ f - all the sets Q with 
the same n loi and kf give rise to the same portfolio for 77. Hence, if we call T the set of 
all these sets, we just have to sum up over all these different contributions: 

III. 2 Second order corrections 

We start considering the easiest situation, i.e. the case of a time independent perturbation 
Hf. the prices are constant in time. Hence the integrals in formula (I3.12p can be easily 
computed and the result is the following: 



where 



?(t) ^//r..y(())//r,/;(())c':,-,;, /■.(/). (3.21) 

T 

Ejr - EjT Q I E Tf - E Tq Ejr f - Ejr 



Recalling definition (13 . 13[) . we rewrite equation (I3.2ip as c^ (t) = < ^t^H-i^t >< 
(fFjHjtpjTQ > £r,ro,r s if) which explicitly shows that up to this order in our perturbation 
expansion transitions between states which differ, e.g., for 2 shares are allowed: it is 
enough that some intermediate state ip^r differs for (plus) one share from (pjr Q and for 
(minus) one share from (fjr . 

If the P a (£)'s depend on time the situation is a bit more complicated but not very 
different. Going back to Example 2 above, and considering then a simple (but not trivial) 
situation in which the prices of the shares really change, we can perform the computation 
and we find 

cP f (t) = {-i) 2 Y,{ h FfAto)h^ {to)Jo{F,J r o,Ff,t 1 ) + 
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+hjr f ^(t 1 )(h J r^ (t )Io(J 7 , JF ; ti)/i(^>, T\ t 2 ) + hjr^it^J^J 7 , F , T f \ t 2 ))+ 



(3.22) 



where we have introduced the functions 



Iji^g-t) := / e 



i{E T -Eg) 



1 



,i{Er-Eg)t _ i(E T -E g )tj 



) 



and 



dt'. 



Needless to say, this last integral could be explicitly computed but we will not show here 
the explicit result, since it will not be used. 

The same comments as above about the possibility of having a non zero transition 
probability can be repeated also for equation (13.221) : it is enough that the time-depending 
perturbation connect ip^ to <p^ via some intermediate state yjj in a single time sub- 
interval in order to permit a transition. If this never happens in [0, t], then the transition 
probability is zero. We can see the problem from a different point of view: if some 
transition takes place in the interval [0, £[, there must be another state, (fjn, different 
from (pr f , such that the transition probability Pj? _>j7> (t) is non zero. 

III. 3 Feynman graphs 

Following [11] we now try to connect the analytic expression of a given approximation 
of CpAt) with some kind of Feynman graph in such a way that the higher orders could 
be easily written considering a certain set of rules which we will obviously call Feynman 
rules. 

The starting point is given by the expressions (13.1 ip and (I3.12p for c^ f (t) and Cj?(t), 
which is convenient to rewrite in the following form: 



c«(i) = -i f 
Jo 



i 



(3.23) 
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and 

4]{t) = {-if V f dt 2 T dt ie tE ^ < ip^Hfafa > x 
T Jo Jo 

x e^* 1 < cp^Hjih)^ > e~ tE ^ (3.24) 

A graphical way to describe Cj?(t) is given in the figure below: at t — to the state 
of the system is tpj? , which evolves freely (and therefore e - *^' 1 ^^ appears) until the 
interaction occurs, at t — t\. After the interaction the system is moved to the state 
(fFf, which evolves freely (and therefore e~ lE:F f tl (fjr f appears, and the different sign in 
( I3.23P is due to the anti-linearity of the scalar product in the first variable.). The free 
evolutions are the upward arrows, while the interaction between the initial and the final 
states, < ifjr^H^t^ip^ >, is described by an horizontal wavy line. Obviously, since the 
interaction may occur at any time between and t, we have to integrate on all these 
possible ti's and multiply the result for —%. 




Figure 1: graphical expression for c)p (t) 

In a similar way we can construct the Feynman graph for c^\t), c^\t) and so on. For 
example (t) can be deduced by a graph like the one in Figure 2, where two interactions 
occur, the first at t = t\ and the second at t — t^: 
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Figure 2: graphical expression for c)p (t) 

Because of the double interaction we have to integrate twice the result, since t\ £ (0, t-z) 
and £2 £ (0,t). For the same reason we have to sum over all the possible intermediate 
states, jr. The free time evolution for the various free fields also appear, as well as a 
(— i) 2 . Following these same rules we could also give at least a formal expression for the 
other coefficients, as cp(t), Cj?(t) and so on: the third order correction c^ (t) contains, 
for instance, a double sum on the intermediate states, allowing in this way a transition 
from a state with, say, n° a shares to a state with n[ a = n° a + 3 shares, a triple time 
integral and a factor (— i) 3 . 

IV Conclusions 

We have shown how quantum statistical dynamics can be adopted to construct and an- 
alyze simplified models of closed stock markets, where no derivatives are considered. In 
particular we have shown that both the Schrodinger and the Heisenberg pictures can be 
successfully used in the perturbative analysis of the time evolution of the market: how- 
ever, the approximations considered in the Heisenberg picture are not completely under 
control because of the many assumptions adopted as it happens using the Schrodinger 
wave function of the market and looking for transition probabilities. 

We have also shown that the Feynman graphs technique can be adopted for the pertur- 
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bative analysis of our market, and some simple rules to write down the integral analytic 
expression for the transition probabilities have been deduced. 

Of course a more detailed analysis of the model should be performed, in particular 
looking for those adjustments which can make more realistic the market we have de- 
scribed so far. These should include for instance source and sink effects to mimikate non 
conservation of the number of shares, short terms exchanges, financial derivatives and so 
on. Incidentally, we believe that the hamiltonian H in (12. 5p and the related differential 
equations in (12 .7p . could be interesting by themselves, hopefully in the description of some 
(realistic?) many-body system. We hope to consider this aspect in the next future. 
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